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Abstract. We study the zero-temperature properties of the Kondo lattice model 
within the dynamical mean-field theory. As impurity solver we use the numerical 
renormalization group. We present results for the paramagnetic case showing the 
anticipated heavy Fermion physics, including direct evidence for the appearance 
of a large Fermi surface for antiferromagnetic exchange interaction. Allowing for 
the formation of a Neel state, we observe at finite doping an antiferromagnetic 
metal below a critical exchange interaction, which shows a crossover from a local- 
moment antiferromagnet with a small Fermi surface for weak exchange coupling 
to a heavy-fermion antiferromagnet with a large Fermi surface for increasing 
exchange. 

Including lattice degrees of freedom via an additional Holstein term we observe 
a significant suppression of the Kondo effect, leading to strongly reduced low- 
energy scale. For too large electron-phonon coupling we find a complete collaps 
of the heavy Fermi liquid and the formation of polarons. 
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1. Introduction 

Heavy Fermion systems based on 4f or 5f intermetallics are paradigms for electronic 
correlations in solid state physics. The low-temperature physics of these compounds is 
strongly influenced by the local moment on the f shell, subject to an antiferromagnetic 
exchange to the conduction electrons. The resulting physical properties are in many 
cases again Fermi liquid like, however with extremely enhanced Landau parameters, 
in particular an effective mass up to three orders of magnitude larger than the one 
found in conventional metals p!Q| l27 l [28]. This large effective mass is the reason why 
these systems are referred to as Heavy-Fermion materials (HF). Moreover, in addition 
to these extreme Fermi liquid properties, the HF also show various phase transitions 
and the thermodynamics show that these transitions actually occur within the heavy 
Fermi liquid [27 . Finally, the appearance of superconductivity in a system with 
initially well-defined magnetic moments is a rather unconventional feature, and close 
investigation revealed early on that the nature of the ordered state may be rather 
unconventional This observation has been substantiated by the development 
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over the past 15 years which showed that a larger number of these HF systems exhibit 
rather peculiar quantum-phase transitions, partially identified as the driving force 
behind the superconducting transitions [28l [29] . 

Other materials which show a coupling between itinerant quasi-particles and 
localized spins are certain transition metal oxides [13], magnetic semiconductors or 
semi-metals in the series of the rare earth monopnictides and monochalcogenides 
[2TJ[26], and diluted magnetic semiconductors such as Gai-ajMn^^As [l4l[T8]. Here, 
the coupling between local spin and conduction electrons usually mediated through 
Hund's exchange and thus typically is ferromagnetic. 

A theoretical description of HF compounds is conventionally based on the Kondo- 
lattice model (KLM) 

The operators c"^^ denote annihilation (creation) operators of itinerant quasi-particles 
with dispersion e^, Si is the operator for the conduction states' spin density at 
lattice site Ri and Si describes a spin of magnitude S localized at site Ri. The 
interaction between the spin of the conduction states and the localized spin is modeled 
as conventional isotropic exchange interaction J. 

The dilute version of the KLM ([T]), the so-called single-impurity Kondo model 
(SIKM) where there exists only one additional spin at site Ri = 0, is well understood 
and shows for antiferromagnetic coupling J < the Kondo effect [12j, which precisely 
leads to the phenomena observed in the Fermi liquid phase of HF systems, viz a 
strongly enhanced mass. There are nowadays several computational tools to treat 
the SIKM, for instance continuous-time Monte-Carlo [20 or Wilson's numerical 
renormalization group [4- 

In theoretical treatments one usually ignores the lattice degrees of freedom. On 
the other hand, all the above mentioned materials have a rather strong electron-phonon 
coupling [T1[T3] and one can expect that the charge physics driven by phonons somehow 
competes with the spin physics due to the exchange interaction with the localized 
spin. Moreover, without exchange coupling, phonons will lead to conventional s-wave 
superconductivity. Thus the investigation of the interplay between a coupling to a 
spin and the lattice degrees of freedom is highly interesting. 

This paper is intended to give a summary of the physical properties of the 
KLM as seen by the dymanical mean-field theory. This necessesarily excludes 
nonlocal phenomena like unconventional superconductivity, but allows for the study of 
antiferromagnetism and whether it is always accompanied by a breakdown of the large 
Fermi surface. Inclusion of phonons eventually leads to superconductivity [3 , which 
however is of the standard local 5- wave type. A detailed account of an investigation of 
the interplay of phonon-mediated superconductivity and HF physics will be presented 
elsewhere. 

The paper is organized as follows. In the next section we discuss the model and 
the approximation used to solve it. The case without phonons, i.e. the conventional 
heavy-fermion physics both in the paramagnetic and the antiferromagnetically ordered 
state is the subject of section [3| The effect of phonons on the low-energy properties 
will be discussed in section [4| A short summary and outlook will close the paper. 
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2. Model and Method 

The KLM Hmiltonian ([T]) will be again the basic model. Except for one dimension, 
no analytical solution exists, and even conventional numerical tools such as Quantum 
Monte-Carlo (QMC) become rather cumbersome due to a severe sign problem away 
from particle-hole symmetry. Thus, a reliable approximate method is needed. If 
one is not interested in the properties too close to a phase transition or in the 
rather complicated, non-local ordering phenomena, a suitable tool is the dynamical 
mean- field theory (DMFT) [9l. Here, the lattice is mapped onto an effective single- 
impurity problem, which can be then solved using standard techniques. Here, we use 
the numerical renormalization group (NRG) approach 01 [30]. One of its apparent 
advantages is the possibility to access small energy scales without problem and cover 
the whole range from T = to finite temperatures of the order of the bare energy 
scales. Furthermore, it also allows to include phonons to a certain extent, namely an 
Einstein mode coupled through a Holstein term to the charge degrees of freedom. As 
is suggestive from effects like Kondo volume collapse [ll [13] , such a term could be 
rather important. We will thus work with a Hamiltonian (for a detailed introduction 
and further references see [24 J 

H = Hklm +u;oJ2blb,+gJ2 (clc^a " l) {4 + ) (2) 

i ia 

where uq is the frequency of an appropriate optical mode and g a measure of the 
electron-phonon coupling. How such an additional local coupling can be treated within 
NRG, is described in detail in [4 . To obtain reasonably accurate spectra also for higher 
energies, we use the broadening strategy introduced by [8]. 

As is well-known, one major effect of the Holstein model is to introduce an effective 
attractive interaction to the electronic subsystem: If the phonon frequency uq and the 
electron-phonon coupling g become large, keeping g'^ /cjOq constant, the phonons can 
be integrated out, yielding an attractive local Coulomb interaction /7eff = —'^g'^/^o- 
Without explicit exchange interaction J, one will then obtain a Hubbard model with 
attractive [/, which shows charge-density and superconducting ordering phenomena 
[21 [7]. From the point of view of Kondo physics, the negative U will lead to a Kondo- 
like behavior in the charge sector, strongly competing with the spin Kondo effect 
introduced by J. We thus can expect interesting physics to occur when both couplings 
are present. 

3. KLM without phonons 

Let us begin with a comprehensive summary of the basics of heavy- fermion physics. We 
will use a 2D square lattice with nearest-neighbor hopping for the conduction states. 
Note that the DMFT is rather insensitive to the dimensionality, and we chose the 2D 
lattice to facilitate visualization of the results. Calculations were done with an NRG 
discretization parameter A = 2, between 1000 and 5000 states were kept per NRG 
step and, where applicable, 50 bosons kept initially. These values were systematically 
changed for selected calculations to ensure that the results are independent of these 
numerical parameters. 
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3.1. Paramagnetic metal 

Let us start with a comparison of the properties at finite J, but with ^ = 0. 
We can distinguish two cases, namely an antiferromagnetic exchange J < and 
a ferromagnetic J > 0. The resulting density of states (DOS) for T = and 




Figure 1. DOS at T = for the KLM ([TJ with \ J\ = 0.25I¥, where W denotes 
the bandwidth of the conduction electrons. The filhng of the band is ric = 0.8. 
The dotted curve displays the DOS at J = as reference. 

I J|/VK = 0.25 at a filling ric = 0.8 of the conduction band is shown in Fig. [l] The 
quantity W denotes the bandwidth of the conduction band and will serve as energy 
scale hereafter. There are notable differences between the two cases J < and J > 0. 
The DOS for J > looks very much like the DOS of the bare conduction band (dotted 
curve in Fig. [T]), although it is somewhat broadened. For J < 0, however, the DOS 
is strongly modified |19J, showing a pseudo gap close to the Fermi energy u = 0. 
This latter feature is a fingerprint of heavy Fermion physics, resulting from a picture 
of hybridized bands [10 . This interpretation becomes even more apparent when one 
looks at the spectral function along the standard k directions in the first Brillouin zone 
of the 2D square lattice in Fig. |2j left part. Dark color means low intensity, bright 
color high. Included as full black line is the bare bandstructure of the 2D nearest- 
neighbor tight-binding band. For J > we basically see the band structure of the 2D 
nearest neighbor tight-binding band. There is a moderate broadening, which actually 
is to be expected, because for J > the local spin effectively acts as a potential 
scatterer for the band states such a situation the DMFT is equivalent to a 

CPA calculation, which yields a constant broadening. For J < (right panel), on the 
other hand, there is a fiattening of the band structure close to the Fermi energy and 
a structure similar to a hybridization gap opens. The fiat portion of the lower band 
corresponds to a large effective mass of the quasi-particles. Note that we have a rather 
sharp structure at the Fermi energy, i.e. one can indeed talk about quasi-particles here. 
Another remarkable difference between the cases J > and J < is observed 
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Figure 2. Spectral functions (left) and momentum distribution (right) for 
I J| = 0.25VK, other parameters as in Fig. ^ The thick line in the lower left 
corner of the contour plot for J < denotes the position of the jump in the 
momentum distribution. 

when one looks at the momentum distribution function n{k) displayed in the right part 
of Fig. |2j As already noted for the spectral functions, the result for J > resembles 
the Fermi function, slightly smeared out by incoherent scattering from the spins. In 
any case, the Fermi surface is located at the kp of a non-interacting 2D tight-binding 
band with a filling Uc = 0.8. On the other hand, the momentum distribution for 
J < does not show any distinct features at this particular value of k. Instead, one 
notes a small jump in n(^) at a vector outside the square marking the Fermi surface 
of a half-filled system. A closer inspection shows that this k vector corresponds to a 
Fermi surface of a system with ric = 1.8, i.e. the system shows a "large" Fermi surface 
with the spin degrees of freedom contributing to the quasi-particles now. Note that 
the height of the jump in n{k) is directly related to the inverse effective mass of the 
quasi-particles. 

The results for J < in Fig. [2] represent the essence of HF physics, namely the 
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generation of heavy quasi-particles, represented by flat bands with a structure known 
from hybridized bands and a large Fermi surface. 

3.2. Antiferromagnetic ordering 

As is well known, the Kondo lattice model shows a variety of magnetically ordered 
states [6l [TTJ [151 [16] . Within DMFT, an overview was presented in [5[[22]. The major 
findings at T = are: At half filling (Kondo insulator regime) one finds a critical 
Jc < 0, with no magnetic phases present for J < Jc and antiferromagnetism for 
J > Jc The staggered magnetization of local spins and band electrons are opposite for 
J < 0, leading to an effectively reduced total moment especially close to the "quantum 
critical" point. Note that this effect is enhanced when one includes interactions in the 
band electron system (see e.g. [22]). 

Away from half filling there appear, at strong enough doping, ferromagnetic 
phases in addition [22^ 25 . For J < Jc, no further magnetic phase was observed in the 
vicinity of half filling. For J > Jc, one can however again stabilize antiferromagnetic 
phases. For example, for a filling Uc = 0.9 of the conduction band, such an ordered 
state is found for J < —W/A. To obtain a reasonable convergence and a stable 
solution, the DMFT calulations must be done with Broyden mixing [31]. Quite 
interestingly, we are not able to find stable and reasonably converged solutions for 
-0.125 < J/W < -0.1. The results for J/W = -0.05 and J/W = -0.15 are shown 




Figure 3. Total DOS for the A sublattice at T = for J/W = -0.05 (full 
line) and J/W = —0.15 (dashed line). The table to the left shows the staggered 
magnetizations. 

in Fig. [3] together with a table of the values for the staggered magetization of the local 
spin, 777^, the conduction band, m^, and their sum, m^. A rather interesting question in 
this connection is whether the ordered state corresponds to a "local moment" regime, 
i.e. where the local spins are effectively decoupled from the band states and one has a 
small Fermi surface, or if it is a "heavy Fermion" magnet with a large Fermi surface. 
The most interesting case of course is when both appear as function of J and there 
might be a phase transition associated with the change in Fermi surface topology. 
Obviously, inspection of the DOS in Fig. [3] alone is not sufficient to identify a possible 
transition. 
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Such information can however be obtained by inspecting the spectral function. 
The results are shown in Fig. [ij where the ^-vectors are now restricted to the first 
magnetic Brillouin zone (see for example [23 ). Included in the figure is as continuous 
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Figure 4. Spectral functions for k from the magnetic Brillouin zone for ric = 0.9. 
Left: Full energy range. Right: Magnified view around the chemical potential. 

black line the dispersion calculated with the Hartree term of the self-energy alone [17] . 
This curve thus represents the proper approximation for the situation of ordered local 
moments polarizing the Fermi sea. The left part shows the overall structures and the 
right part a magnified view of the region around the Fermi energy. 

Quite apparently, the concept of a small Fermi surface, i.e. ordered local spins 
polarizing the band states, describes the situation very well for J/W = —0.05. The 
situation changes completely for J/W = —0.15. First, the staggered magnetization 
has dropped to ^ 0.12 now, with ml ^ 0.2 and ^ 0.08 (see Tab. |3|. Thus, 
while the polarization of the band states has not changed significantly, the one for the 
local spins has dropped by more than 50%. Further increasing \J\ does not lead to any 
new features, the polarizations of both local moments and band electrons smoothly 
drop to zero for J/W ^ —W/A. Second, the spectral function has become much 
more HF like with a flat band around F (c.f. Fig. [2|. Therefore, the scenario of a 
decoupled Fermi sea and local moments ordering does surely not apply any more, in 
fact the whole system looks much more band like, with a large Fermi surfacejf] This 
scenario is also supported by the fact, that for increasing coupling \J\ the system more 
strongly shows tendencies to rather form spin density waves than Neel order. These 

t Note that "large" in the magnetic Brillouin zone means centered around F, as the M point has 
been mapped back to the zone center. 
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observations strongly hint towards a phase transition from a local moment like phase 
at small \J\ to a phase where magnetic ordering appears in the heavy quasiparticles at 
larger \J\. However, a decisive answer what type of phase transition this is at present 
cannot be given. 

4. Effect of phonons 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 

0)/W 

Figure 5. DOS at T = for the KLM ^ with J = -O.bW and p = 0, ^ = 

and g = 0.4VK. The filhng of the band is ric = 0.8. The inset shows the dependence 

of the effective mass on g. 

An extended account of the effect of Einstein phonons on HF physics is given in 
[24] , where we concentrated on the periodic Anderson model. Here, we show further 
results for the KLM, which supplement the findings reported there, in particular the 
behavior as the Kondo effect collapses due to polaron formation. We restrict the 
discussion here to J < and fix the phonon frequency to uq = 0.5VK. Calculations 
were done for J = — 0.51^ at a filling Uc = 0.8. The values for ujq and J are apparently 
rather large. However, smaller J and ujq do not change the qualitative picture, but 
make it much harder to visualize the structures. 

The results are summarized in Fig. |5]for ^ = 0, ^ = 0.21^ and g = 0.41^. The 
first thing to note is that the phonons lead to a reduction of the width of the pseudo 
gap close to the Fermi energy and also a reduction of the overall bandwidth. In 
addition there occur new structures at higher energy with increasing electron-phonon 
coupling, which are related to the formation of polarons. The reduction of the overall 
bandwidth is expected and can be interpreted as an increase of the effective mass 
of the bare conduction states due to the coupling to the phonons. The reduction of 
the width of the pseudo gap, on the other hand, signals a likewise reduction of the 
low-energy scale generated by the Kondo effect. Both effects lead to an effective mass 
as function of g as depicted in the inset to Fig. [5] Note that m* initially depends only 
weakly on g. However it diverges very strongly as g ^ ujq. 
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Figure 6. Spectral functions for J = —0.5W, other parameters as in Fig.js] High 
intensity is represented by red, low by violet. 

These feature becomes more apparent by inspecting the spectral functions in Fig. 
[6j Both effects, the overah reduction of the bare bandwidth and the increased HF 
mass, are clearly visible. Moreover, with increasing g one finds roughly ^-independent 
structures representing the polaronic modes. Further increasing ^, we observe a rather 
sharp crossover around ^ ^ cjq to a completely incoherent behavior. The results of a 




co/W 



Figure 7. DOS (left panel) respectively spectral function (right panel) for 
g = 0.5W = OUQ. Other parameters as in Fig.js] 

calculation for g = ujq are shown in Fig. [7j Note that there is no Kondo feature left 
either in the DOS or in the spectral function, and all structures are rather broad. We 
also would like to mention that for g > ujq it becomes increasingly hard to stabilize a 
given occupation n ^ 0, 1, 2 of the conduction band. This indicates that the system is 
close to a Peierls instability, i.e. the formation of a charge density wave together with 
a lattice distortion. 

5. Summary 

We have presented a summary of properties of the Kondo lattice model within 
dynamical mean-field theory at T = using the numerical renormalization group as 
impurity solver. We have extended the Kondo lattice model by including an Einstein 
mode coupled to the electrons via a Holstein term. The importance of such modes for 
HF materials can be deduced from strong effects such as the Kondo volume collapse 
observed in Ce. 
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Without the phonons, we find for antiferromagnetic coupling the expected heavy- 
fermion behavior, with hybridized bands appearing in the spectrum and a large Fermi 
surface. For ferromagnetic coupling on the other hand, the bare band structure is only 
weakly modified due to incoherent scattering from the local degrees of freedom. 

Allowing for a magnetically ordered state, we are able to stabilize a Neel order 
for J > Jc, where Jc < depends on the filling of the conduction band. For the case 
Tie = 0.9 discussed here, we find Jc/W ^ —0.25. Quite interestingly, there appears 
to be a qualitative difference between the ground-state for| J| and \J\ \Jc\. In 
the former case, the results can well be interpreted with a system consisting of local 
moments which order antiferromagnetically via an RKKY-like exchange. In this case, 
we have a small Fermi surface and the band states are weakly polarized by the presence 
of the local moments. For the latter case we find a rather different behavior. First, 
although the polarization of the local moments is considerably smaller than for small 
J, the band polarization has actually increased pointing towards a much stronger 
entanglement between local spins and conduction electrons. Second, the spectral 
functions again form fiat bands around points of the magnetic Brillouin zone which 
map to the heavy fermion bands without Neel order. A more detailed investigation of 
the crossover respectively transitions could not yet be accomplished due to convergence 
problems of the DMFT in the interesting region. 

Adding phonons, we find a general narrowing of the bare band, which also leads 
to a reduction of the Kondo scale. Eventually, when the coupling becomes of the 
order of the phonon frequency, the electrons tend to localize and form polarons with 
the phonons. At that point, the effective mass diverges and the electronic spectrum 
becomes incoherent. As a side observation we note that in this region one sees a 
tendency of the system to form a charge density wave. 

The results for T = presented here strongly motivate further investigations, 
in particular at finite T, searching for the critical temperature and in particular the 
interplay between HF and local moment physics with respect to magnetism, charge 
odering and superconductivity. 
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